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Abstract 

A collision of orbifold planes in eleven dimensions has been proposed as an expla- 
nation of the hot big bang [1, 2, 10]. When the two planes are close to each other, the 
winding membranes become the lightest modes of the theory, and can be effectively 
described in terms of fundamental strings in a ten dimensional background. Near 
the brane collision, the eleven-dimensional metric is an Euclidean space times a 1+1- 
dimensional Milne universe. However, one may expect small perturbations to lead into 
a more general Kasner background. In this paper we extend the previous classical 
analysis of winding membranes to Kasner backgrounds, and using the Hamiltonian 
equations, solve for the wave function of loops with circular symmetry. The evolution 
across the singularity is regular, and explained in terms of the excitement of higher 
oscillation modes. We also show there is finite particle production and unitarity is 
preserved. 
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1 Introduction 



The initial singularity problem remains an open question in Cosmology and any model of 
the early Universe requires a resolution of this paradigm. We know general relativity breaks 
down close to it but there is hope that a theory of quantum gravity can resolve the singu- 
larity. Recently, interest has turned to the particular case of bouncing models, where the 
question of how information propagates across the big crunch/big bang transition has not 
been completely solved. For example, in the ekpyrotic/cyclic model this transition is as- 
sumed smooth with controlled particle production [1, 2]. In order to prove or disprove such 
a statement, most people have used effective field theories, constructed from string theory 
or other extensions of general relativity. There is plenty of literature on this approach to 
tackle the singularity problem (see for example [3]), however, most of these effective theories 
break down near the singularity, and perhaps one should be considering a more fundamen- 
tal description beyond general relativity, such as string/M-theory. One such approach is to 
directly investigate the string equations of motion in a singular background [4], with special 
attention being given to the Milne universe [5, 6]. Some authors have used this to argue that 
this particular singularity can not be resolved [6], however, the results are not conclusive and 
in fact evidence from a dual description such as investigated in [7] seems to contradict the 
result. Furthermore, there is evidence that in a big crunch/big bang transition in asymp- 
totically Anti de Sitter spacetime, the Conformal Field Theory (CFT) description leads to 
a well-defined evolution of fields across the singularity [8, 9]. In [10], the authors proposed a 
novel approach to explain the Milne singularity using eleven dimensional membranes, which 
is a natural setup for the cyclic universe, as also discussed recently in [11]. In this paper, 
we generalize this M-theory setup to more general backgrounds, corresponding to the ho- 
mogeneous and anisotropic Kasner metrics. Such a background could well result once we 
include the effects of small perturbations in the background isotropic metric. Moreover, we 
also make progress in the quantum evolution of such membranes across the singularity. 

In the M-theory model of [10], the singularity is described by two orbifold planes that 
collide as the eleventh dimension, which separates them, disappears. The study focuses on 
the evolution of membranes stretching from one orbifold plane to the other, but in particular, 
considers winding membranes, which correspond to the lowest Kaluza Klein modes in ten 
dimensions. As argued in [10], the winding membranes represent the lightest modes and de- 
couple from the bulk (heavy Kaluza Klein) modes when the eleventh dimension is sufficiently 
small. From the ten dimensional point of view, these winding membranes are described by 
perturbative string theory, hence they include perturbative gravity. The classical evolution 
has been studied in the case of the Milne universe [13], and some progress has been made 
to understand the quantum theory of such modes, either by taking semiclassical approxima- 
tions, such as the instanton calculations of [10] or by prescriptions to linearize the classical 
equations of motion, as in [13]. Furthermore, as shown in [14] the classical evolution does 
not acquire finite-width a'-corrections either far away from the singularity or very close to it. 
Therefore, on either side of the singularity there are two semiclassical regimes connected by 
a phase where quantum corrections are important. In the present analysis the story repeats, 
but in this case we quantize the action for certain membranes - corresponding to circular 
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strings from the ten dimensional theory - and show how particle production remains finite 
even though higher oscillation modes are excited. 

In the ekpyrotic/cyclic universe, the spacetime is well described by the Milne universe 
near the orbifold collision [1, 2]. However, it is well known that as the singularity is ap- 
proached, any small perturbation to the Milne universe can lead to a Kasner solution (see 
for example [12]). Therefore, we believe it is important to show how these winding mem- 
brane modes evolve across Kasner metrics. The metric, g^, for an eleven dimensional Kasner 
space-time is given by 

10 

ds 2 u = -dt 2 + \doA 2n {dx 1 ) 2 , (1) 
i=i 

where 8q is a dimensionful positive constant, which for the eleventh direction represents the 
rapidity at which the two orbifold planes collide. The usual Kasner conditions hold in eleven 
dimensions: 

10 10 

£rf = l = X> ( 2 ) 

i=i i=i 

We have chosen the singularity to be at t = and we have glued the manifolds before and 
after the singularity using the absolute-value function. In general, t = is a curvature 
singularity, and only for the particular case of pi — 1 for a given coordinate, does the 
solution become a direct product of a 9d flat space-time and the Milne universe, with t — 
a coordinate singularity, which simply represents the fact that we have made a bad choice of 
coordinates in flat space. However, if the spatial coordinate of the Milne metric is compact, 
then the singularity is a conical singularity. This is the case of the cyclic universe where 
the big crunch/big bang transition is modeled by an orbifold collision, where the eleventh 
dimension is compact with a Z 2 symmetry. The orbifold structure is not essential for the 
present discussion, because our results only rely on a very small compact eleventh dimension, 
so we will forget about this discrete symmetry. In other words, because we are concentrating 
only on the Bosonic sector of the theory, where all the string models have the same field 
content, our results apply to either Heterotic or II A limit of M-theory. 

When the eleventh dimension, x 10 , is small enough we can use a ten dimensional descrip- 
tion based on the Kaluza-Klein reduction 

ds 2 n = e- 2 ^ 3 ds 2 10 + e 4 ^ 3 d(x 10 ) 2 , (3) 

where the dilaton is given by = In \ 6 t\ (with p 10 = p^), and the ten dimensional metric 
reduces to 

ds 2 w = a{tfds 2 conf = a{tf (^-dt 2 + \0 t\ 2p '(dx^ , a(t) = \9 t\^ 2 , (4) 

where we have assumed p^ > 0, implying the eleventh dimensions disappears as t — > 0. As 
x 10 we can think of these winding membranes as fundamental strings on the orbifold 
planes feeling the metric (4). Alternatively, the winding membranes can be thought of as 
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strings with a time- dependent tension living on the metric ds 2 con ^ as will become evident 
later. Since the string coupling is e^, for really small times — close to the singularity — the 
strings hardly interact and one can take the free string action as a good description. There- 
fore, we will focus our attention on the propagation of free strings on the ten dimensional 
dilaton-Kasner background (4). 

The paper is organized as follows: in Section 2, we write down the different actions for 
membrane excitations in eleven dimensions, and in the following section we solve the classical 
equations of motion governing only winding membranes with cylindrical symmetry. Section 
3 is devoted to the quantum description of circular loops using a Hamiltonian approach, 
before we finally conclude in Section 4. 



2 Winding membranes 

Our starting point is a Polyakov type of action for a bosonic membrane of tension /i 2 in 
eleven dimensions 

Spd = -y J d 3 aC po i (5) 

= -y J d 3 ay/^ {l^d^d^g^ - l) , (6) 

where x M are fields representing the position of the membrane in a target space with metric 
g^y. The worldvolume spanned by the coordinates a a has a metric 7 Q/ g, and the variation 
of this action with respect to 7 a/ 3 yields the constraint 7 Q( g = d a x^dpx l 'g^, which can be 
substituted back into the action to obtain the Nambu-Goto action, 

Sng = -A*2 J d 3 a^J-Det(d a x^d^g^). (7) 

The first action is more convenient to analyze the quantum behavior whereas the second is 
more useful to describe the classical evolution, as we will show below. As explained in [10], 
the Hamiltonian can be constructed from the action (6), leading to the constraints 

U = n^g^ + f j%Det{d&x' i d $ x v g lll ,) =0, V & = n^x" = 0, (8) 

where 7r M = ^§0- are the canonical conjugate momenta to x^, x^ = |^ and the hatted 
indices run over the spatial dimensions of the membrane's worldvolume. Therefore, the 
most general Hamiltonian is 

H = Jd 2 a (^H + A & V^j , (9) 

where the two functions A and A a represent the gauge freedom of the membrane's metric 
diffeomorphisms. We consider a partial gauge where the momentum is always orthogonal to 
the membrane, which is equivalent to choosing A a = 0, and will use the remaining gauge 
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freedom to simplify the equations of motion and obtain either classical or quantum solutions. 
A winding membrane is obtained by demanding its coordinates x^ are independent of one of 
the spatial membrane worldvolume coordinates (say a 2 ), except for the eleventh dimension 
which should be proportional to a 2 . We choose x 10 = a 2 (where a 2 runs from to 1), so 
that after integrating with respect to a 2 in (7) we get an overall factor of |6*ot| P0 in front of 
the effective string action. 



3 Classical evolution 

To describe the classical evolution of a winding membrane in a Kasner background we use 
the t = t gauge in the Nambu-Goto type of action (7). Then the action reduces to 



dadT\0 o tf* 



1 - J2 \dot\ 2p '(i 1 ) 2 J Yl IW^idaX*) 2 , (10) 



1=1 



where r = er° and a = a 1 . For simplicity we will assume = 1 during the calculations and 
then restore a general 9q at the end. The equations of motion, in units of /i2 — 1, read 



X* = — — — — - — — — , 7h = \t?^ + ^d, ( ' " X ' 



e V e 

„i\2 



+ 2p i )\t\ 2{p " +2pi) ^ L , A = < + \t\ 2(p * +2pi) (d a a^) 2 , 

t Si 



where the string energy density 7r = e is given by 

e 2 = X>|- 2 ^ 2 - (12) 

i=i 

Using the last expression we can rewrite the differential equation for ^ in the following way 

q ( ^ \ = 2(^ + 2p t ) 

which will be useful later. From (11), divergent solutions arise when at least one of the 
Kasner exponents in the lOd space-time is negative enough to lead to a divergent term in 
the energy density e at t = 0, as previously shown by Tolley [15]. On the other hand, regular 
solutions across t — are obtained if all pi > To avoid divergences we will assume 

Pi > -p<t>/2 for all i. (14) 

The divergent cases correspond to situations where, before the singularity, one spatial di- 
mension expands faster than the contraction of the lOd conformal factor a(t), as appreciated 
in (4). We are more interested in situations which are small perturbations away from the 
Milne universe, but still close to it. 
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To construct a perturbative solution around the singularity, one can expand the equations 
of motion (11) in terms of the string tension, as was done in [13]. Formally, one introduces 
a parameter A in place of the tension (i.e. fi — > fi\) and solves iteratively the equations of 
motion as a series in A. At the end, one sets A = 1. As shown in [13], the only equation where 

A appears is -hi = \\t\ 2( - p 't' +p ^ d a (j^f-^j , whose solution to zeroth order in A is 7Tj = 7Tj(0), where 

7Tj(0) = 7Tj(0,(T) is the loop momentum at t — 0. Assuming (14) holds, one can integrate 
equation (13) and insert the solution of e and 7Tj in the ij-equation, to obtain the zeroth 
order solution 

f* 7r(Q)\t \~ 2pi 
x l (t)~x l (0) + dt tKm + C(A), (15) 

Jto JEj \t\' 2P3 h 2 (°) + (^'(o)) 2 \t\ 2p ^] 

where x l (0) = x l (0, a) is the string shape at the singularity. The last integral is finite, hence 
the solution is regular at t — 0. In general, this integral has to be done numerically, but there 
are specific string geometries or configurations where the solution can be found analytically; 
this is the case of a circular string. 



3.1 Circular loops 

For the rest of the paper we will focus on the circular loop, which as explained in [13] is the 
classical analogue of the dilaton field. The simplification rests on the fact that the only dy- 
namical coordinate is the radius of the circle. Furthermore, to preserve the circular symmetry, 
the two Kasner exponents of the plane where the loop oscillates should be equal. Without 
loss of generality, we assume the circular loop oscillates in the xy plane and has a center of 
mass velocity v in the z direction, with the ansatz x % = (R(t) cos(cr), R(t) sin(cr), vt, 0, 0), 
and the Kasner exponents in these directions are p = pi = p 2 and p 3 = p z . Under these 
assumptions the equations of motion (11) simplify to 



7T 



U\-2p z +p 

V = — = , 7T Z = 0, 

e 



R=^tl, *R = -\t\*»**-, 

e I 

I = { P4> + 2p)\t\ 2 ^ +2p ^, e 2 = 7T 2 + \t\ 2 ^K 2 z + \t\ 2 ^ +2 ^R 2 , (16) 

where ir 2 R = irf + -k\ and e = t p e. Again, we can rewrite the differential equation for q using 
the last constraint, namely 

a ( e ~ 2 ^ - 2 (^ + 2 p) f:r 2 uvfp-p,)^) (17) 

Moreover, in the case of a circular loop it is not hard to find another constraint by combining 
the different equations in (16), given by 

(y 2 + R 2 )\t\ 2p + 4 |t|- 2 fr*+P> = 1. (18) 
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Notice that the speed (squared) of any point in the loop is V 2 = V 2 + Vr = (v 2 + R 2 )\t\ 2p , 
which is unaffected by the contraction or expansion of the plane of oscillation. After a careful 
analysis of the second term in the last constraint, one can be convinced that every point 
in the string reaches the speed of light (V 2 — >■ 1) as t — >■ 0, if the inequalities (14) hold. 
As a result of this, the solutions are not time invariant and the outgoing mode is different 
from the incoming one. Quantum mechanically this time asymmetry is the origin of particle 
production and excitation of higher order oscillation modes, as we will show later. This effect 
is enhanced when the center of mass momentum n z vanishes. Furthermore, in the case when 
the bound (14) is saturated for the p Kasner exponent, the speed of the loop will be finite 
and generically smaller than the speed of light; hence there will not be a time asymmetry in 
the solution and no particle production or higher oscillation modes will be expected across 
the singularity. This should be expected, since the effective metric on the xy plane (see 
equation (4)) neither contracts or expands when the bound (14) is saturated for p. 

Although the set of equations (16) can not be solved analytically everywhere they can be 
solved approximately in different regions and these results can then be compared to the full 
numerical solutions. The solutions evolve similarly to those in [13], so we refer the reader 
to this previous work for details. However, we would like to stress a few general points, 
especially when more general Kasner exponents are considered and not only the Milne case, 
as it was done in [13]. 

Far away from the singularity, the string does not feel the contraction or expansion of the 
universe, and therefore, it oscillates as if it lived in the ds 2 con f metric of equation (4), which 
reduces to flat spacetime for the Milne universe. Following the notation of [13], at a time 
to the winding membrane in eleven dimensions can be effectively described by perturbative 
string theory in ten dimensions. Furthermore, the string coupling = | 1 3p0//2 tends to zero 
as the singularity is approached, hence the free string action becomes more accurate closer 
the orbifold collision. By definition, to corresponds to the time where the string coupling 
hits unity, namely 

to = do 1 - (19) 

The string tension is therefore /ii = |6> to| p< ^ 2 = 1 and its length is l s ~ ^ l — 1. In terms 
of a quantum analysis we would expect this regime to be well described by a semiclassical 
solution, which can be obtained using the WKB approximation. However, classically we 
start with a string configuration in the metric ds 2 con j at time t = —t and evolve it towards 
the singularity. The solution crosses the singularity and after it has reached a large enough 
positive time (comparable with t = +t ), we can trust the description of strings living on 
the metric ds 2 on j- again. We aim to compare both states, the ingoing and the outgoing states 
at t — ±t respectively, to determine whether there was particle production or excitation 
of higher vibrational states. Quantum mechanically, it corresponds to calculating a mini 
S-matrix, defined by the evolution of the quantum ingoing states at t — —t to the outgoing 
ones at t — +t . 

Starting from a large negative time in this adiabatic vacuum (t = —to), the string evolves 
into the singularity increasing its size, according to the contraction of the conformal factor 
of the universe. By a simple rescaling of the coordinates it is possible to show that the 
loop scales as a power of l/a(t) ~ t~ Pr t>/ 2 . Once the size of the loop is comparable to the 
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averaged conformal Hubble radius (~ the "stringy" quantum corrections become 

really important and the evolution can no longer be described by a semiclassical analysis. 
Remember the string coupling goes to zero as t — > 0, hence only the a'-corrections become 
important as the solution approaches the singularity. As shown in [14], these a'-corrections 
modify the semiclassical evolution, but in a finite way. Finally, close to the singularity, 
there is another semiclassical phase in which the modes freeze, stop oscillating and cross the 
singularity. In detail, the string "breaks" into string bits which evolve independently of each 
other. In other words, the spatial gradients that tie the string together become negligible and 
the evolution only depends on time. This phenomenon is a consequence of the ultra-local 
behavior that one expects near a cosmological singularity (see for example [12]). At t — 0, 
however, the string receives an energy kick, because it has to travel at the speed of light, 
which either increases or decreases the amplitude of the outgoing mode, leading to classical 
gain or loss of energy. Quantum mechanically this translates into particle production, as we 
will discuss in the next Section. 

Different Kasner exponents can dramatically change the outcome because they play an 
important role close to the singularity. In order to simplify the discussion, we will only 
consider the n z = case, but the analysis can be easily generalized for a non-zero center 
of mass momentum. Therefore, from the equations of motion (16) one may consider the 
exponents p and p^ as two free parameters of the model, however, the lid constraints (2) 
relate these parameters, together with the exponents of the orthogonal directions. To simplify 
the argument, one may think of the orthogonal exponents as being split into two sets, the 
first being all the same (m of them) and the other set being zero (6 — m of them), i.e. Pi = p 
(i — 3, m + 3) and pj = (j — m + 4, 9). Then, after solving for p^ in terms of p and 
m using the constraints (2), one obtains 

V<i> = [(1 - 2p) ± y/m(m + 4p - 2(m + 3)p 2 )}/(m + 1). (20) 

To get real exponents, p should lie between (1 — a/1 + m(m + 3)/2)/(m + 3) < p < (1 + 
y/l + m(m + 3)/2)/(m + 3). There is therefore a wider range of allowed values for p as m 
increases, but also divergent solutions are more likely to arise since p<f, may be smaller than 
— 2p, and especially for the negative branch of the p^ solution above (see Figure 1). 

To have a feeling of how sensitive the solutions are to the Kasner exponents let us vary 
p around the Milne solution (p = and p^ — l). Figure (2) shows a slight variation around 
the Milne universe, and even though the general behavior remains similar, the amplitude 
and periodicity of the outgoing modes strongly depends on the precise value of the Kasner 
exponents. Figure (3) shows a more dramatic change near the singularity when the Kasner 
exponents are taken to be far away from the Milne universe. If p ~ —p^/2, the energy density 
has a very mild dependence on \t\, so the modes do not feel the contraction or expansion 
of the universe. Therefore, there is no classical energy production, and as we will see later, 
there is also a mild particle production in the quantum theory, which consistently tends to 
zero for p = —p^/2. The opposite case, when p is positive and relatively large, the string 
feels effectively a larger contraction/expansion of the scale factor, which produces a big effect 
around t = 0, and thus a greater quantum production of particles. 
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Figure 1: Kasner exponent p<f> as a function of the Kasner exponents in the xy plane (p) and 
the m equal-valued Kasner exponents in the orthogonal directions. The blue solid straight 
line is for m — 0, the solid line ellipsis is for m = 1 and the blue dotted one for m = 6. 
Divergent solutions which cannot be followed across the singularity correspond to values of 
p<f) below the red dashed line which only exist for m > 1. 



As mentioned before, we can describe the evolution of these strings as an expansion in 
the string tension around t = 0. Formally, we introduce a parameter A where the string 
tension is and then truncate the expressions to the desired order in A. Finally we set A = 1. 
This was done above explicitly for a general string configuration to zeroth order in the string 
coupling, and in the case of circular loop (with tc z = 0) the integral (15) can be solved 
analytically to give: 



Sign(7T^) 
R = R ° + t (l-p)\t\ P 2Fl 



1 — p 1 1 — p R t 



2 



2(p^ + 2p)'2' 2(p + 2p)'7r° 



+ 0(X), (21) 



where R and ir° R are the values of R and n R at t = 0, and 2 Fi is the Gauss hypergeometric 
function, defined as the series 

2 Y 1 [aMc ] z] = Y^^f^z\ (22) 

with (w) n = w(w + l)...(w + n — 1) for any complex number w. The hypergeometric function 
is well behaved for all values of the Kasner exponents that satisfy the bound (14). Solution 
(21) reduces to that found in [13] for p^ = 1 (and p = 0). 



4 Quantum description 

We now turn our attention to the more complicated problem of quantization. If one tries 
to naively quantize the classical equations (11), all sorts of problems arise, because of the 
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Figure 2: The radial coordinate R and the energy density e evolving in time t for different 
Kasner exponents around the Milne solution (jp^ = 1 and p — 0). The evolution is regular 
across t — 0, when the energy density takes its lowest value. The solutions are not time 
reversal and strongly depend on the Kasner exponents. We have assumed 7i z = 0, and the 
initial conditions are R(—t ) = 1 and R(—t ) = 0, with t = 20. 



square root present in the action. However, one can take a different approach, by considering 
the Polyakov type action and using the Wheeler-de Witt formalism, namely 

H$> = 0, (23) 
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Figure 3: The radial coordinate R evolving in time t for different Kasner exponents far away 
from the Milne solution (p^ = 1 and p — 0). For positive p the kick around the singularity 
is stronger and the outgoing solution is drastically modified. On the other hand, negative 
values for p, close to saturating the bound (14), lead to relatively little modification of the 
outgoing solution, resulting in less particle production, as will be explained in the quantum 
model. We have assumed 7r 2 = 0, and the initial conditions are R(—t ) = 1 and R(—t ) = 0, 
with t = 20. 



where \I> is the wave function of the string. However, it is hard to proceed from here, 
because the Hamiltonian includes a term proportional to d a x° = d a t, so it is difficult to 
synchronize the different points on the string and to talk about a common time for the 
string. Fortunately, if we restrict ourselves to circular strings where the only degree of free- 
dom is the radius of the loop, R, then the Hamiltonian simplifies enough for the problem 
to be tackled. In order to preserve the circular symmetry in time, we also need the Kas- 
ner exponents of the plane where the loop oscillates to be equal, as we assumed in the 
previous section. Moreover, a circular loop can have an initial center of mass momentum 
perpendicular to the plane of oscillation. We take the circular winding membrane ansatz 
x^ = (t(r, a 1 ), R(t) cos(ct 1 ), R(t) sin(cr 1 ), vr, 0, 0, cr 2 ), where without loss of generality we 
allow the string to oscillate in the xy plane, with pi = p 2 = p, and a constant center of 
mass velocity v in the z direction with Kasner exponent p 3 — p z . Notice that the circular 
symmetry forces d c t = 0, which implies t is a function of r only. Then the Hamiltonian (9), 
reduces to 

H = J d 2 a^ [-(tto) 2 + \9 t\~V R + \9 t\- 2 ^7r 2 z + \9 t\ 2 ^R 2 ] , (24) 
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where, as mentioned before, we have chosen the partial gauge A a = 0. To simplify this 
expression we can perform a canonical transformation to a new time coordinate given by 
t ~ t 1+p , up to a constant factor. Finally, because the Hamiltonian density does not depend 
on the spatial worldvolume variables, we can choose the extra gauge freedom to fix A = 
2t 2p /(J d 2 a), simplifying the Hamiltonian to 

H = -(n ) 2 + tt 2 r + \6 t\ 2 ^-^7r 2 z + \9 t\ 2 ^ +2 ^R 2 , (25) 

where ft = t p 7i is the canonical momentum of t, and the new Kasner exponents are 

P4> = P<t>l( l +P)i P = P/(1+P), (26) 

where we assume p ^ — 1 (which is consistent with (14)). In fact, it is possible to rescale t 
and R (preserving the Poisson brackets) in such a way that 9$ only appears in the tt z term. 
Then, one can redefine the center of mass momentum to absorb the information of both, 
the orbifold rapidity 9 and the loop's center of mass momentum n z , into a single combined 
parameter, 7r , in the following way 

n e = 9 { * Ap+s)/(s+2) n z , (27) 

where 

s = p <p + 2p, Ap = p- p z . (28) 

Notice that making 9 larger increases rr c if Ap > —s/2, which is true within some open set 
of the allowed parameter space of the Kasner exponents. Evidently, this open set includes 
the Milne universe (p = p z = and p^ — 1). So in order to understand the physics of the 
cyclic model with respect to the orbifold rapidity, in what follows, we will consider Kasner 
combinations which obey 2Ap > —s/2. 

However, solutions of the wave function using the Hamiltonian (25) are messy when the 
term with ir c is present, thus we will focus our attention to the limit 7r c — > first and then 
present the results for a non-vanishing ir c . One can neglect 7r c in two different limits: either 
when the orbifold rapidity is very small compared to the loop's momentum, or when the 
center of mass velocity is small enough on its own 1 . 

4.1 Zero center of mass momentum 

When dropping the center of mass momentum term, some observables, such as the amount 
of particle production, do not depend explicitly on the orbifold rapidity, which is given by 
6*o. This statement may sound like a contradiction because we would expect more particles 
being produced for larger collision rapidities, but if one looks in more detail, the result is 
consistent with QFT in curved space. The basic idea of particle production is to measure 
the "difference" between two vacuum states, which in our case correspond to one in the far 
past, before the singularity, and another in the far future. However, a change in 9 not only 

1 Note that 9 < 1 is needed to start with loops inside the Hubble horizon during the contracting phase, 
as shown in [13], so that in this case tt z ~ implies tt c ~ (where we have also assumed 2Ap > —s/2). 
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changes the orbifold rapidity but also the two vacua we are comparing, and for the case of 
7r c — > the changes are such that their effects cancel each other. 

Therefore, knowing that we can absorb 8q by a constant rescaling of the variables in the 
case of negligible n c , we will set 9 = 1 and then restore it in the final expressions. Thus, for 
the present discussion we consider the simpler Hamiltonian 

H = -(n ) 2 + n 2 R + \t\ 2s R 2 , (29) 

A canonical quantization (7r — > idi and n R — > iOr) implies that the Wheeler-de Witt 
equation (23) reduces to 

[dt ~d 2 R + Itl 2 ^ 2 ]^, R) = 0, (30) 

where we have dropped the tilde over the time variable for simplicity, but the reader should 
remember that we are actually still referring to t defined earlier. Eqn. (30) corresponds to 
the Klein-Gordon equation with the potential of a harmonic oscillator with a time-dependent 
frequency u(t) ~ \t\ s . There is a parallel line of thought to understand this result: as was 
mentioned earlier, an alternative description of a string evolving across the metric (4) is to 
think of a string oscillating in the space-time ds 2 on j- with a tension that varies with respect to 
time (i.e. fii = /i2^o^| P0 )- From this point of view, a natural wave equation for the string is 
that of a harmonic oscillator in the metric ds 2 con j with a time-dependent mass, which after a 
change of variables can be recast into a time-dependent frequency problem (see for example 
[16]). Furthermore, when solving string equations in singular backgrounds, we generally find 
harmonic oscillator equations with a time-dependent frequency (see for example [17]). 



4.1.1 Asymptotic solution 

The first step to construct a solution to equation (30) is to understand the asymptotic 
behavior, where the adiabatic or WKB approximation provides a good approximation to the 
true solution. Naively we would expect this to be in the adiabatic regime when \t\ is large 
enough, however, as we will see later, the story is a bit more complicated and the adiabatic 
regime holds only for large T time (T ~ £( 2 + s )/ 2 ) 2 . Nevertheless, we can still get a feeling of 
the asymptotic solution if we set the frequency to be constant, namely 

[d?-d 2 R + t 2 s R 2 }^(t,R) = 0, (31) 

where t » 1. The wave function that solves this simplified equation is related to the usual 
harmonic oscillator form, and is given by 

oo 

ty(t, R) = Y^ H n(t S o 2R ) ex PHo#7 2 ) [Aexp(iE n t) + B ex$(-iE n t)\ , (32) 

n=0 

where E\ = (2n + l)tg = (2n + l)w(t) and H n ( x) is the Hermite polynomial of degree n. 
Notice that the energy levels E n are actually the square root of the usual harmonic oscillator 

2 Remember the effective (2+f )-metric where the string oscillates is given by equation (4), namely 

d 4d = \t\ p H-dt 2 + \t\ 2 Pdx 2 + \t\ 2 Pdy 2 ) = \t\ (P * +2p) \-dP + dx 2 + dy 2 ), where t - t*> +1 and p and p given 
by (26). Therefore, the "10-dimensional" time T is simply given by T = J dT - / t s ' 2 dt, with s = p + 2p. 
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levels. To understand this, one can think of the classical string analysis previously done, 
where in fact, the Hamiltonian e (see equation (16)) is the square root of the harmonic 
oscillator for constant time. 

The asymptotic solution (32), provides a hint as to the best ansatz we can adopt to 
obtain the general solution to (30). Using the harmonic oscillator as a basis and replacing 
t by t, we consider an ansatz of the form: 

oo 

*(t,R) = J2 A n(t)H n (\t\ s/2 R)exp(-\t\ s R 2 /2), (33) 

n=0 

where the A n 's are considered to be functions of time and determine the evolution of the 
harmonic oscillator states from the incoming vacuum at negative times to the outgoing modes 
far after the singularity. We can then use the orthogonality of the Hermite polynomials 

J H n (x)H m (x)e- x2/2 dx = ^{Tn\)5 mn (34) 

and the relations H' n (x) = 2nH n -i(x) and H n+ i(x) = 2xH n (x) — 2nH n ^\(x) to decompose 
equation (30) into an infinite system of coupled ODEs for the A n 's, which are given by the 
following recursive equation 

= A n -^A n +[4s-s 2 (l + 2n + 2n 2 ) + lQ(l + 2n)\t\ 2+s ]-^ (35) 

+(n + l)(ra + 2)jA n+2 - (2 + s){n + l)(n + 2)^A n+2 - -^i n _ 2 

2 2 

+ (2 + s )-^2 A n-2 + (n + l)(n + 2)(n + 3)(n + 4)^p4»+4 + g^A>-4- 

One should note that these equations are not regular at t — 0, so one cannot find solutions 
which interpolate between negative times and positives ones. This implies nothing else than 
using the harmonic oscillator basis is not a good approximation around the singularity. 
However, as we will see in the next section the evolution across the singularity is simpler 
than one could possibly have expected. 
Consider the zeroth mode A Q equation, 

Ao ~ Y t A ° + ( 4s - s " + 16 N 2+S ) A) + jA 2 - ^|^2 + ^-A 4 = 0. (36) 

If we assume the higher modes are negligible, A n = for n > (which should be the case if 
the incoming state is the vacuum and will be justified later), then Eqn. (36) can be solved 
in closed form for either positive or negative times, resulting in 

A = \t\( 2+s » 4 (9 ) s/4 [ Cl ^ (1) (2A^ /2 |t|( 2+ ^ 2 ) + c 2 Kl 2) (2k6f \t\^' 2 )] , (37) 

where (x) and Kf 1 (x) are the Hankel functions of the first and second kind respectively, 
ci and c 2 are integration constants, k = 1/(2 + s), / = A; a/1 + s 2 /2, and we have restored 
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6*o in the expression. Using the asymptotic expansion of the Hankel functions Ka' 2 \x) = 
^7r/(2x) exp (±i(x — (2a ± we obtain an asymptotic solution for the ground state 

given in terms of positive and negative frequency components, namely 

V(t, R) = exp(-\6 t\ s R 2 /2) [d exp(iE \t\) + C 2 exp(-iE \t\)] , (38) 

where C\ = \J\c\ exp(— i(2a + l)7r/4), C 2 = a/T c 2 exp (i {2a — l)ir/4) and now the frequency 
is time- dependent and given by Eo(t) = ^uo(t) = 

To justify dropping the higher order modes (A n with n > 0) in equation (36) for large 
times one can re-write the recursive differential equation in terms of the T time ( T = 
_|_^s/2^(2+s)/2^ an( j ftien see which are the leading terms for large T. Equation (35), in 
T-time and with 9 — 1, reads 

= ^ + (l + 2n)A n + O^y (39) 

and has a large T solution, 



A n (t) = exp(iv^TTT) = exp(iE n t), E n (t) = |t| s/2 . (40) 

The picture is then the following: asymptotically (large T), all the harmonic oscillator modes 
decouple from each other and as one approaches the singularity they start interacting. If 
one starts with the ground state as the incoming vacuum, then higher order modes become 
excited in order to resolve the singularity, and one ends up with a tower of states as the 
outgoing state in the far future, when again the interaction of harmonic oscillator modes 
stops and we can use this basis to describe the resulting state. Furthermore, since the in- 
equation (35) is a "double-step" recursive equation, if one starts with the ground state in 
the far past, then only even n-levels will be excited across the singularity. 

4.1.2 Solution near the singularity 

If we consider the vacuum positive energy mode (38) as the incoming state for t — > —to 
(with t » 1), then we can evolve numerically the solution to the full equation (30). The 
numerical solution behaves regularly everywhere, particularly at t — 0. In order to picture 
the evolution of the wave function across the singularity, we can use the integrated wave 
equation 

/+oo 
V(t,R)dR. (41) 
-oo 

The numerical solution of such an integrated function is shown in Figure 4, and notice in 
particular that around t — the solution is simply given by a straight line, as we shall 
shortly describe. This result provides a simple explanation of how the modes evolve in terms 
of the absolute value function. The fact that the wave function (41) is well approximated 
by ty(t) = at + b around t — 0, can be understood as follows: assume ^ behaves like 
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Figure 4: Evolution of the integrated wave function ^(t) in time for the Kasner exponent 
combination p^ = 0.793521 and p = 0.4. The positive frequency mode in (38), with C\ = 1 
and C*2 = 0, at t = —20 was chosen. Different initial times do not change the behavior 
but the overall scale. Around t = the wave function evolution follows a straight line, and 
after the singularity, higher order modes contribute to the overall wave function, leading to 
a non-periodic and complicated structure. 



ty(t,R) = e~^ SR2 / 2 f(t) near the singularity 3 . After substituting it into the Hamiltonian 
equation (23), and assuming a simple polynomial function for f(t) of the form f(t) = at m + c, 
then the only solution with a well-defined limit as t — > has c = and m = s/2, (2 + s)/2. 
Therefore, the solution near t = 0, which is valid on both sides of the singularity, is 



which after integrating over R implies \l/(t) = at + b. If we think of the harmonic oscillator 
basis, the functions A n should all behave as [a n t + b n )\t\ s ^ 2 for small times on either side of 
the singularity. However, it doesn't mean we can simply solve for the individual coefficients 
a n and b n , because from (33) and (41) we require knowledge of all the modes that cross 
the singularity, thus no individual mode can be evolved across using this description for the 
integrated wave function around the singularity. Fortunately, from the numerical solution, we 
can always calculate the A n coefficients using the orthogonally property (34) of the Hermite 
polynomials. Therefore, we can show the agreement between the numerical solution and 
the expected behavior ((a„t + b n )\t\ s ^ 2 ). For example, Figure 5 shows the case of the zeroth 
mode A before the singularity. 

3 If one plots the wave function near t = as a function of R, one finds a good fit using the Gaussian 
profile exp(— \t\R 2 /2), which justifies our assumption. 




(42) 
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Figure 5: Using the orthogonality property (34) one can extract the coefficients of the 
harmonic oscillator modes A n from the numerical solution of the wave function (see equation 
(33)). Here, we plot Re(A ) for p^ = 1 and p = 0. The dashed line represents the analytic 
solution (37) which follows the numerical solution (solid line) up to the place where the 
t = behavior, given by solution (42) becomes more accurate (dotted line). We normalize 
the wave function so Aq ~ 1 initially. 



On the other hand, after t = the .R- integrated wave equation ^(t) no longer looks 
periodic (see Figure 4). However, all of this complicated structure can be decomposed as a 
tower of higher harmonic oscillator modes, which were excited during the transition through 
t = and converge as we increase the mode number (see Figure 6). The choice of Kasner 
exponents has an effect on the tower of excited states, as will be shown later when calculating 
the particle production. 

Furthermore, if we follow the whole evolution of the vacuum state, the adiabatic regime 
is well approximated by the analytical solution (37), even as we approach the singularity. 
However, such an agreement inevitably breaks down very close to the singularity, as the 
expansion in terms of uncoupled harmonic oscillator states breaks down. That is, all higher 
modes begin to get excited in order to solve through the singularity, hence, we can not 
neglect the A2 or A\ terms in the Aq equation near t = 0. However, notice that these higher 
order terms never dominate, as can be seen in Figure 6. As the adiabatic approach fails, the 
above mentioned polynomial behavior near t = takes over and the solution can be followed 
all the way to and across the singularity. Some time after t = 0, the story repeats and 
we can use the adiabatic description again, leading to a decomposition of higher oscillation 
state, which can be described in terms of Bogoliubov transformations, as we will show in the 
following section. 
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Figure 6: Using the orthogonality property (34) we plot the A n 's for different mode numbers 
n. One can see how higher modes get excited after t = 0, and the amplitude of higher states 
decrease with n as expected. We use p^ = 0.793521 and p = 0.4. 



4.1.3 Particle production 



One has to be careful when constructing the Hilbert space, since a C°° and globally hyper- 
bolic manifold is needed. In our case the string metric is not well defined at t = 0, therefore 
we can split the space of solutions into two sections: for t < and t > 0. Based on the reg- 
ularity of the classical and quantum solutions around t = it is natural to assume unitarity 
at t — 0, which allows us to use the Bogoliubov transformations to relate the two sections of 
the Hilbert space and calculate particle production. Our numerical results will act to justify 
this assumption. Moreover, a natural inner product in this Hilbert space is defined by [18] 



(*i,* 2 ) = 



dR 



dt 



(43) 



One can expand the wave function ^(R, t) in terms of two basis sets, one for asymptoti- 
cally negative times and one for the corresponding positive times, namely 



*(R,t) = ^^^ + (4^^)* t^-OO 
n 



(44) 



where a' n and a n are the creation and annihilation operators associated with the harmonic 
oscillator expansion, and which obey the usual commutation relationships (i.e. [a n ,ajj = 
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5 mn , etc.). We would expect the negative time vacuum (cin n ^\0)i n = 0) to be different from 



that of positive times (ai?^ \0) out = 0), leading to particle production. The functions U n ' 
and Un ' ae> are the positive frequency modes 



(in) 



= N n exp(+iE n t)exp(-\t\ s R 2 /2)H n (\t\ s/2 R) t < 0, 
jjiout) = N n exp(-iE n t)exp(-\t\ s R 2 /2)H n (\t\ s/2 R) t > 0, 

with the normalization factor M n given by 

M n = (2" +1 n! v / (2n+l)7r)- 1 / 2 , 



(45) 



(46) 



in such a way that (\I/ n , \l/ m ) = S nm . 

Actually, since we are interested in having the vacuum as our incoming state, the wave 
function for t — > — oo takes the form 

* = Mo [(at } + (af V) cos(£ |t|) - i (af } - (af V) sin(S |t|)] exp(-|t| s J R 2 /2), (47) 
where Ao = (4^174- Now, we can write explicitly the wave function for t — > +00 as 



* = E-A/; (^ + (0^)008^ 

n L 

= [ (4^ + (^"V) COs(E |t| + ¥>„) 



/2 T?\-t°R 2 /2 



H n {t s ' 2 R)e 



-1 (af } -(afV) P n sin(E |t| + ^) 



/ 2 f?^-* SR2 /2 



H n (t s/2 R)e 



(48) 



where £> n and </? n (£> n and </? n ) are the amplitude and phase — with respect to t — - 
of the outgoing mode n after sending an incoming cosine (sine) piece of the incoming wave 
function (47). Since the two lines in the previous equation are equal, the coefficients should 
be related in the following way 



a n ^ ~ a n0 a ^ + Aio( a O ^> 



(49) 



where the Bogoliubov coefficients are defined as 
a n0 = ^{V n e^+V n e^), 



(3no = ^{V n e^-V n e 1 ^). (50) 



Particle production at a given energy level can then be understood in terms of the expectation 
value of the particle number operator N n = (a^^On 1 ^ over the out-vacuum |0) OMt , which 
simplifies to the following 



(0|iV n |0> ont = f3 n0 f3* n0 = T- 2 n\V2^TT {V 2 n + V 2 n - V n V n cos( Vn - <p n )) . (51) 
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Figure 7: Bogoliubov coefficients a^ and (3^ as functions of the mode number n for three 
Kasner exponent combinations: p = 0.4, p = 0.0 (Milne) and p = —0.37, with p^ obtained 
using (20) for m — 7. In the last case, the effective (3+l)-metric where the string oscillates is 
almost flat, thus the particle production \(3\ 2 is negligible. A fit shows that both Bogoliubov 
coefficients decay as C\ exp(— C^n 3 / 4 ), where C\ and Ci are positive constants that depend 
on the Kasner exponents only. 



Figure 7 shows the particle production for different Kasner exponents. A fit of this plot 
shows that the particle production decays exponentially with mode number, n, satisfying 
a fit of the form /3q exp(— Cin 3 ^), where /?oo and C\ depend on the particular choice of 
Kasner exponents. Of particular note is that for large n where the energy per mode, E n 
satisfies E\ ~ n, (see equation (40)), we then obtain an exponential decay of the form 
/3oo exp(— CiEn 2 ), a result that agrees with the semiclassical instanton calculation in [10] for 
the Milne universe, but remains to be compared for other Kasner exponents. 

4.2 Non-zero center of mass momentum 

If one starts with a non- vanishing center of mass momentum, it z , at t = —to, and if the 
orbifold rapidity 6$ is not negligible with respect to 7r 2 , then one cannot drop its contribution 
in the Hamiltonian (25). The Wheeler-de Witt equation associated with the Hamiltonian 
(25) can be solved numerically, again yielding a regular solution, in particular at t — 0, 
which is very similar to the tc c = case shown in Fig. 4. Furthermore, to describe the 
t > behavior, we can once again use the harmonic oscillator basis as an ansatz, but with 
a frequency which now depends on both time and 7r c ; details are given in Appendix A. Here 
we only summarize our findings. Asymptotically, the wave function can again be described 
in terms of decoupled harmonic oscillator modes, while near the singularity the straight-line 
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behavior (after integrating over R) remains as a characteristic feature. Thus, we can send in 
the vacuum state and read off the harmonic oscillator components in the outgoing solution, 
to find the Bogoliubov coefficients for such a transition. 

From the definition in equation (27), we know that the orbifold rapidity should play a role 
in the physical observables, via the combination ix Q = #Q 2Ap+s - ) / (s+2 - > 7r 2 . Thus, if as mentioned 
above we assume Ap > —s/2, increasing the orbifold rapidity is equivalent to increasing the 
center of mass momentum. Now the question is: how do the Bogoliubov coefficients change 
due to 7r c ? On the one hand, we know that classically the presence of a singularity forces the 
loop to increase its velocity until it reaches the speed of light at t — 0. As we have already 
seen, this speeding up is achieved by exciting higher order oscillation modes. Therefore, 
one would naively expect that a larger n c should help the loop to reach the speed of light 
at t — without exciting higher order modes. On the other hand, one should expect that 
increasing the orbifold rapidity makes the orbifold collision more violent, thus increasing 
particle production. These two effects compete against each other and indeed, we see both 
of them at different scales. For small n c , particle production is suppressed exponentially 
as exp(— C\ix\ — c 2 n 3 / 4 ), where C\ > and c 2 > depend on the Kasner exponents, and 
c 2 has a mild dependence on n c . However, for larger n c particle production grows as a 
powerlaw of 7r c , while exponentially decreases with n. This powerlaw growth is in agreement 
with the instanton calculation of [10], where it was found that {0\N \0) out oc ^" 1)/3 for 
a (d + l)-dimensional spacetime. In agreement, we find for the case d — 4; three spatial 
dimensions corresponding to the effective space in which the loop oscillates plus the M- 
theory dimension. Figure Fig. 8 shows how particle production decays with n, and Fig. 9 
depicts how the particle number of the ground level, {0\N \0) out , changes with n c . At first 
sight, it may be worrying to think that the particle production diverges for large 7r c , but 
actually one should have in mind that the orbifold rapidity should be small, in order to have 
modes created well inside the comoving Hubble horizon where the adiabatic regime can be 
trusted [10, 13]. Furthermore, the transverse velocity of any loop cannot be the speed of 
light, otherwise it could not oscillate in the xy plane. 

For a gas of loops, the average velocity should be close to l/\/2, as shown in [19] for flat 
space 4 . We can then imagine a single loop with such a velocity at t — —t s (see equation 
(19)), and calculate the corresponding value of n c using the first equation in (16) and (27). 
Assuming the loop is initially static (R(—t s ) = 0) and of unit size (R(—t s ) = 1), one obtains 



v6, 



(2Ap+2)/(s+2) 







- VT^ ■ (52) 

If 6*0 is small and the velocity dispersion around the mean value is not large, then the particle 
production should be small. For example, for 9 Q ~ 1/5 and a velocity dispersion of around 
±0.25 ([20]), the corresponding values of ir c lie on the range e G (0.3, 1.93) (at t = —t s = —5), 
which maps to the exponentially decaying region of Fig. 9. 

Finally, an important check of the quantum model is to verify unitarity, which implies 
that the canonical commutation relationships should be preserved over time. This translates 



4 We can assume a flat background initially, since for a small 6q the size of the loops is small compare to 
the Hubble radius . 



21 



0.8 
B § 0.6 - 
0.4 
0.2 




i 1 r 

(Milne) 

X 7l c =0.0 

■ 7l c = 1.0 

+ 7l c =7.0 

♦ ji_= 14.0 



0.4 - 



0.1 



(Milne) 

7l c = 0.0 
Jt c = 1.0 
7l c = 7.0 
7C C = 14.0 



Figure 8: Bogoliubov coefficients and /3^ as functions of the mode number n for four 
different initial momenta, 7r z , measured at t = —t Q = —10. The Milne universe (p = 0, 
= 1) is assumed for all cases. Fits show that both coefficients decay as c 2 exp(— cin 3//4 ), 
where C\ has a mild dependence on 7r c , but c 2 strongly depends on 7r c as shown in Figure 9. 
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Figure 9: Bogoliubov coefficient /3q for different initial initial momentum 7r c chosen at t s = 5, 
and assuming the Milne universe (jp — 0, p$ — 1). For small center of mass momentum the 
loop needs less excited modes to reach the speed of light at t — 0, but a large 7r c is equivalent 
to a large orbifold rapidity, which increases the particle production due the more violent 
brane-collision. The fitting curve has a form c 2 exp(— C\ir1) + csn^, where the q's are positive 
constants. 
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into the following property for the Bogoliubov coefficients 



nO- 



(53) 



One would need to sum over all states to get unity, however, one can see that using the first 
ten excited modes the number is close to unity (in Figure (10)). 
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Figure 10: Plot of £ n = E?=o(«?o 
three Kasner metrics of Fig. 7, and for two non-zero values of 7r c at t s 
preserved if J2 n ~ > 1 as n — > oo. 



5. Unitarity is 



5 Conclusions 

The nature of the big bang is one of the biggest problems facing cosmology. No one fully 
understands the physics that occurred in this crucial period. One of the more interesting 
recent proposals for the origin of the big bang, is through the collision and re-emergence of 
two orbifold planes in eleven dimensions [1, 2, 10]. As the two planes approach each other, 
the light states of the theory consist of winding M2-branes, which are described in terms of 
fundamental strings in a ten dimensional background. Near the brane collision region, the 
full eleven-dimensional metric considered is that of Euclidean space times a compactified 
1+1-dimensional Milne universe. In [13], two of us considered the the classical evolution of 
winding membranes in such a background, showing that they suffered no blue-shift as the 
M theory dimension collapses, and their equations of motion remained regular across the 
transition from big crunch to big bang. However, one may expect there to really be small 
perturbations to the background metric, leading to a more general Kasner background for 
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the light membranes to evolve in. If this is the case, an obvious question is what happens 
to collapsing light M2-branes as they pass through the singularity? In this paper we have 
gone beyond to original classical analysis of winding membranes in [13] to include general 
Kasner backgrounds. By considering the corresponding Hamiltonian equations, we have 
been able to solve for the wave function of loops with circular symmetry and demonstrate 
the sensitivity of the solutions to the values of the Kasner exponents p around the Milne 
solution (p — and p^ — 1). As is evident from Figures 2 and 3, although the general 
behavior remains similar, in particular the loop solutions remain perfectly finite, there is 
clear evidence that the amplitude and periodicity of the outgoing modes depend strongly 
on the precise value of the Kasner exponents. In the regime p ~ —p^/2, the energy density 
has a mild dependence on \t\, and the modes do not feel the contraction or expansion of the 
universe, leading to no classical energy production. On the other hand, for large positive 
values of p, the string does feels the effects of the evolving scale factor, which produces a large 
effect around t — 0, and a greater classical production of energy. To confirm the particle 
interpretation of this classical result, we adopt the Wheeler-de Witt formalism to quantize 
the system of the evolving circular loop in the Kasner background. Circular symmetry of 
the loop is crucial in order to help us solve the system of equations as it allows us to assume 
d a t = 0, implying the time coordinate of the loop is a function of r only. Remarkably we 
are able to solve the Wheeler-de Witt equation both far from and close to the singularity 
t = 0. Asymptotically, the loop is well described through a WKB approximation, given by 
Eqns. (33-38) in the case of very small orbifold rapidity compared to the center of mass 
momentum, or by Eqns. (55)- (57) for the general case. Of particular note though is the 
fact that near the singularity, the solution simplifies to such an extent that the integrated 
wave function (41) is given by a straight line in t through the singularity which allows a 
simple understanding of the solution on either side of t — 0, made evident through Figures 
4 and 5. The complicated non-periodic evolution seen in Figure 4 just after the singularity 
provides evidence that there could well be particle production as a loop evolves through a 
singularity. This is confirmed in our analysis of the particle production, as seen in Figures 7 
and 8 where we plot the particle production number for three separate values of the Kasner 
parameter including the usual Milne case, and three different center of mass momenta. In 
all cases, particle production is exponentially suppressed with the oscillation mode number 
n, and in the case of zero-center of mass momentum tt z it is independent of the orbifold 
rapidity 9q. However, for a non- vanishing center of mass momentum, particle production 
depends on the simple function 7r c = 7r c (9o, 7r z ), defined in (27). Production of particles is 
exponentially suppressed for small n c , but has a powerlaw growth for large 7r c . The two 
effects are expected: the first one corresponds to the fact that a small initial velocity helps 
the loop to reach the speed of light at the singularity without exciting higher order modes. 
In contrast, the second effect corresponds to the fact that a large orbifold rapidity induces a 
more violent big crunch/big bang, which results in more particles been produced. 

As we have seen, as the loop passes through the singular point, higher order modes 
become light and excited, the string reaches the speed of light everywhere along it, and the 
effect is either an increase or decrease in the amplitude of the outgoing mode leading to 
classical gain or loss in energy and the production of particles whilst maintaining Unitarity. 
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There are of course caveats to what we have done in this paper. In reality the Kasner 
metric is not flat away from the singularity and thus it can receive corrections when interact- 
ing with the membrane. It is not clear how this would affect the background and subsequent 
analysis? Assuming we are only slightly away from Milne case, we believe that the analysis 
should be similar to what we have done here. 

It is very encouraging that we have seen how it is possible to have finite particle produc- 
tion in through a singular region in a Kasner background. It is now worth seeing quite how 
much reheating occurs in such a scenario. An obvious, if somewhat difficult calculation is to 
extend this work beyond the case of a simple circular loop. Finally, it will be interesting to 
see how these results can be interpreted from the conformal theory description of [9]. 
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A Non-vanishing tt c solution 

To solve the Wheeler-de Witt equation derived from the Hamiltonian (25) for a non- vanishing 
7r c term, one can use the same ideas as in the n c = case. First and to get a felling of the 
solution, one could study the wave equation derived from (25) for a constant time t = t , 
namely 

H^j = [d 2 t -d 2 R + tf X 2 + P S R 2 ) * = 0, (54) 

where 7r c , Ap and s are defined in (27) and (28). A solution to this equation is given 
again by the harmonic oscillator solution (32), but with an energy E\ = (2n + l)tg + TT^t p . 
Therefore, we can use the same ansatz (33) and the properties of the Hermite polynomials 
to get a recursive differential equation for the function A n (t), which is similar to (35) but 
with extra terms containing n c . Actually, in the T-time variable, the A n equation reads 

= ^ + (1 + 2n)A n + ^±1^^ A n + O (±) , (55) 

which does not have a simple analytical solution for a generic Kasner exponent combination, 
unlike the 7r c = case. In the particular case of the Milne universe, the solution to this 
equation is given in terms of Whittaker functions, which are related to the polylogarithm 
functions. The appearance of these functions suggests a consistent description with the 1 / a'- 
series solution of classical evolution near t — 0, where polylogarithm functions were found 
to second order in the string tension [13]. 

Since we are searching the asymptotic spectrum of states, it is enough to find a series 

/ 2 \ 2Ap— s 

solution of equation (55). The expansion parameter is Q = 7rM ^\T\^ j = 7r;?|t| 2Ap ~ s , 

and the solution becomes more accurate for large |T| because (14) implies — 2Ap + s = 
p^ + 2p z > 0. This series solution can be constructed using the ansatz 

A n = exp{tS n T), (56) 
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where 

£ n = £ n (T) = y/p + aQ + ^Q*. (57) 

The following constants 

P = 2n + 1, a= S + \ , 7 = 2 S - ^ , (58) 

solve equation (55), with an error of order tM 4A p+ s + 2 )/( s + 2 ) ^ ^-(4Ap+s+2)/2^ g Q ^ we can uge 
this approximate solution and calculate the particle production using the same expressions 
as in section 4.1.3. We assume the two vacua are given by 

n 

= ^a^ ut) U^ ut) + {a { ° ut) )\U^ ut) y t^+oo (59) 

n 

where 

= N n exp(+iE n t)exp(-\t\ s R 2 /2)H n {\t\ s/2 R) t < 0, 
jjiout) = N n exp(-iE n t)exp(-\t\ s R 2 /2)H n (\t\ s/2 R) t > 0, (60) 

and E n = j^£ n \t\ s / 2 . The normalization factor N n is given by 

M n =[2 n+l n\^{£' n T + £ n )]- l/ \ (61) 
with £' n ee j^£ n (T). Therefore, by writing the vacuum solution as 

^ — No [(4 m) + (at V) cos(£ |t|) - i (at ] - (af V) sin(£ |*|)] exp(-|t| s J R 2 /2), (62) 

we can read off the outgoing solution in terms of the Bogoliubov coefficients (50), and 
calculate the particle production in the same way we did for the 7r c = case. The results 
are shown in Figure 8 and summarized in Section 4.2. 
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